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We present a unified approach for calculating the properties of shallow donors inside or outside 
heterostructure quantum wells. The method allows us to obtain not only the binding energies of 
all localized states of any symmetry, but also the energy width of the resonant states which may 
appear when a localized state becomes degenerate with the continuous quantum well subbands. 
The approach is non-variational, and we are therefore also able to evaluate the wave functions. 
This is used to calculate the optical absorption spectrum, which is strongly non-isotropic due to 
the selection rules. The results obtained from calculations for Si/Sii-zCe^ quantum wells allow 
us to present the general behavior of the impurity states, as the donor position is varied from the 
center of the well to deep inside the barrier. The influence on the donor ground state from both the 
central-cell effect and the strain arising from the lattice mismatch is carefully considered. 
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I. INTRODUCTION 

Impurity states in heterostructures have been the sub- 
ject of detailed investigations during the last two decades. 
Traditionally, impurities are considered to be purely 
detrimental, by increasing the scattering rates. On the 
other hand, doping is essential to supply enough free car- 
riers into the system. The strive has therefore been to 
remove the doping region from the active region by using 
modulation doping. 

Recently, however, impurities have been placed in the 
active region of Si/SiGe quantum well structures, to ex- 
ploit their presence and properties for novel optical de- 
vices in the far-infrared or terahertz (THz) regioni 2 ^ 
Si and Ge are non-polar materials with low intrinsic 
absorption at THz frequencies. Taking into account 
also the possibilities for integration with existing device 
technology's these systems appear very attractive for op- 
tical applications in the THz-region, which recently has 
received a lot of attention in a variety of fieldsi* A de- 
tailed knowledge of the properties of impurity states in 
Si/SiGe heterostructures is therefore essential. 

Much effort has been spent on calculating the 
impurity energy states in a quantum well (QW). 
The techniques, all within the effective-mass ap- 
proximation, range from the commonly employed 
variational techniques**^**?*****!**** to direct inte- 
gration of the Schrodinger equation** and basis 
expansions pi*** The systems considered have, to the 
best of our knowledge, exclusively been III-V materi- 
als such as GaAs/ALrGai-^As quantum wells, for both 
donors 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 14 and acceptors j 12 i 13 i 15 i 16 and the re- 
sults are found to agree reasonably well with experimen- 
tally obtained valuesi***** 9 ^**** Particular aspects such 
as image charges 9 and the effect of screening of the impu- 
rity potential by the free carriers in the wellSa have also 
been investigated. 

To allow for a more detailed comparison of the cal- 
culated ground state binding energy with experimental 
results, one must also take into account the central-cell 



effect, or the chemical shift. Its physical origin is the ab- 
sence of screening by the valence electrons at distances 
smaller than the outer electronic shells, in which case 
effective-mass theory breaks down. In principle, there- 
fore, a proper treatment of chemical shifts is not possible 
within the traditional effective-mass theoryi** Neverthe- 
less, several authors have considered the central-cell effect 
in the bulk by introducing various short-range pseudo- 
potentials, which are adjusted to give agreement with the 
experimental binding energiesi****i2* This approach has 
also been attempted in the quantum well case^ whereas 
Mueller et al. have chosen to relate the QW central-cell 
shift to the envelope function amplitude!* 7 - 

The donor ground state naturally appears below the 
lowest QW quantization level. As we will see, there are 
however Rydberg-like series of impurity states attached 
to each QW subband, as well as to the three-dimensional 
continuum. The situation is thus very similar to what has 
been observed in strong magnetic fields, where impurity 
states are attached to each Landau leveli** 

It has been noted by several authors that in narrow 
quantum wells, the lowest anti-symmetric impurity state, 
which is bound to the second QW subband, may become 
resonant with the continuum of the first subbandA* Res- 
onant states of this so-called Fano type, i.e. a discrete 
(and hence localized) state degenerate with a contin- 
uum, are well-known from atomic physics 29 and bulk 
semiconductorsi 24 i 30 i 31 ' 32 i 33 i 34 i 35 If symmetry allows it, 
the localized state can couple, or hybridize, with the con- 
tinuous states. The resonant state is then characterized 
by an energy width, which is immediately related to its 
life-time^* It is of crucial importance to determine this 
width in order to consider the effects of the coupled im- 
purity states on the optical and electrical properties of 
modulation-doped quantum wells. 

When the impurity is placed outside the well, in the 
barrier, similar but entirely different resonant states may 
also be formed from the usual shallow donor levels, and 
their widths can be evaluated by the resonant coupling 
methodi* 6 - Attempts have been made to apply an essen- 
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tially equivalent method for the case of donors inside the 
we ^ > n 4 37 J 38 j-j. ^g^g however been shown that the approxi- 
mations made in this approach are rather severe, and the 
widths do not agree with more exact results obtained us- 
ing the same method as will be presented in this paper, 39 
Only when the coupling between the impurity state and 
the quantum well is weak, can such perturbation meth- 
ods be expected to yield accurate results. 

In this work we shall present a non-variational ap- 
proach for calculating the energy levels of shallow donors, 
in which we expand the total Hamiltonian in a complete 
basis. This turns the Schrodinger equation into a matrix 
eigenvalue problem, which is diagonalized to yield all lo- 
calized, hybridized and continuous states. The approach 
has several benefits over the variational method. First of 
all, no assumptions are made regarding the form of the 
impurity wave functions, but instead we obtain the cor- 
rect envelope function, which then can be employed for 
calculating e.g. optical matrix elements. Second, the en- 
tire energy spectrum is treated simultaneously, and hence 
we also obtain a description of the perturbed continuum 
and can observe the formation of the hybridized resonant 
states, and calculate their energy widths. Additionally, 
the difference in the effective mass in the well and bar- 
riers is easily included, as well as an electric field, and 
we are able to consider (but not limited to) anisotropic 
materials such as Si and Ge. It is furthermore possible 
to place the impurities in the barrier instead of in the 
well, and we thus have a unified approach for treating 
impurities in modulation-doped heterostructures. 

We will specifically consider Si/Sii^Ge^ quantum 
wells, which are strained due to the lattice constant mis- 
match. In result, the sixfold degeneracy of the conduc- 
tion band bottom, and the donor states, is partially lifted. 
Moveover, the impurity ground state in bulk Si is known 
to be shifted and split by the central-cell effect 4Mi We 
shall therefore also present an approach for simultane- 
ously taking the strain and the central-cell effect into 
account. 

One of the advantageous properties of a coupled reso- 
nant state is the possibility of populating it by electrically 
pumping the carriers, from the bottom of the subband, 
until they reach the resonance energyiiS 4 ^ 3 . They may 
then be captured into the resonant state, and possibly 
make an optical transition to the ground state. If an 
inverted population could be arranged between the reso- 
nant state and the ground state, one could realize a laser 
based on this process. A particularly appealing point 
of such a device is that since the impurity states are at- 
tached to the QW levels, it is possible to control the intra- 
impurity transition energy by varying the QW parame- 
ters. Experimental evidence of optical transitions involv- 
ing coupled resonant states in quantum wells exists from 
both Raman scattering 4 ^ and absorption spectroscopy 45 
measurements. We therefore also calculate the impurity 
absorption spectrum, for arbitrary polarization. 

The remainder of our paper is outlined as follows. In 
Section^we present the non-variational method for solv- 



ing the Schrodinger equation of a shallow donor in a QW. 
Once the matrix eigenvalue problem is solved, the char- 
acterization of the eigenstates becomes an important is- 
sue, and this is discussed in Section ITTT1 The influence 
of strain and the central-cell effect is considered in Sec- 
tion ^V] and in Section [3 we calculate the quantum well 
optical absorption spectrum in the presence of coupled 
resonant states. The numerical results of the calcula- 
tions for Si/Sii-^Gea; quantum wells are presented and 
discussed in Section I VII followed by a summary in Sec- 
tion ehi 

II. THE BASIS EXPANSION METHOD 

We shall consider the problem of a shallow donor, de- 
scribed by the impurity potential V c (r), in a quantum 
well. Possibly an electric field £ along the QW growth 
direction z is also present. This field may be externally 
applied or built-in due to charge redistribution in the 
structure. In the effective-mass approximation, the total 
Hamiltonian can be written as 

H = H QW + V c (r) + e£z, (1) 

with the quantum well Hamiltonian 46 

h 2 ( d 1 d I , \ . . , , 

#QW = ~ rT"5~~ + ^T V 2D +V(z). 2 

2 \oz m±(z) az mn\\{z) J 

From a principal point of view, our method allows for 
the quantum well potential profile V(z) to vary arbitrar- 
ily with the position z. We will however, for simplicity, 
choose V(z) to be a constant band offset U outside the 
well (which has width a) and zero inside it. The effect of 
strain due to the lattice mismatch between the well and 
barrier regions is assumed to be incorporated in the band 
offset, except for the splitting of the donor ground state. 
This will be considered in detail in Section IIVI along 
with the central-cell potential which is not included in 
the effective- mass Hamiltonian 

We wish to apply our method to anisotropic materials 
such as SiGe alloys. We have therefore split the kinetic 
energy operator into two terms, and let the effective mass 
depend on the direction, with m± and TO|| denoting the 
respective masses perpendicular and parallel to the two- 
dimensional (2D) QW plane. Further note, that both 
the parallel and perpendicular masses are taken to de- 
pend on the coordinate z, due to the presence of the het- 
eroboundaries. In result, the Schrodinger equation with 
the Hamiltonian i?Qw will not be separable. 

The dependence of the parallel mass m\\ on z is usually 
ignored in heterostructure calculations, since the elec- 
trons involved in in-plane transport are strongly confined 
to the QW channel. In the present case, this simpli- 
fication is not appropriate since - especially for narrow 
wells - a large portion of the impurity wave function may 
appear in the barriers, and the binding energy depends 
critically on the exponentially decaying part of the wave 
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function. The discontinuity of the parallel mass can alter- 
natively be taken into account by introducing an effective 
well depth, which depends on the in-plane momentum, 47 
but this is not convenient for our purposes. By diago- 
nalizing IZqw within our basis expansion method (see 
below), we find that a difference in the parallel masses 
leads to non-parabolicity of the QW subbands, but does 
not mix the different QW levels. (This is not to be con- 
fused with the non-parabolicity that arises from coupling 
to the valence band; c.f. the discussion after Eq. ©.) 

The impurity potential is taken as the Coulomb poten- 
tial 



V c (p,z) = - 



(zVp 2 + (* - zo) 2 



(3) 



(in CGS units) as is appropriate for a shallow donor lo- 
cated at (0,0, zq). Here p is the radial vector in the 
QW plane. It has been demonstrated that the discon- 
tinuity in the dielectric constant e is of importance in 
GaAs/Al^Gai-^As systems^ To take this into account 
properly, one could introduce image charges, but to eval- 
uate the matrix elements with the resulting effective im- 
purity potential (c.f. Eqs. (5) and (6) in Ref. 0) would 
lead to rather extensive calculations in our present ap- 
proach. Since the difference in the dielectric constant is 
small between Si and Si/Sii_ x Ge x for small x, we will 
assume that e is a constant, independent of z. 

Our method for solving the Schrodinger equation 
H*S>(r) = E^(r) with the total Hamiltonian (QJ, is based 
on expanding the total wave function ^(r) in a complete 
basis, and diagonalizing the Hamiltonian in this basis. As 
the basis it is natural to use the quantum well eigenstates 

|<?k> = ^ (A'|<?k)=<Kk-k')<W, (4) 

normalized as indicated, where q enumerates the QW 
levels and k is the wave vector for the in-plane motion. 
These states diagonalize 
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where mjj 7 is the effective mass inside the quantum well 
for the direction parallel to the QW plane. 

The Hamiltonian Hq differs from i?qw m that the 
parallel effective mass m\\ does not depend on the z- 
coordinate. Hence, in contrast to the case with -ffq>W) the 
Schrodinger equation with the Hamiltonian Hq is separa- 
ble. To be complete, the basis must include both bound 
and unbound states, and we therefore enclose the sys- 
tem in a box of width C z in the z direction. The box 
is chosen large enough that it has negligible influence on 
the results. Nevertheless, the basis is complete for any 
size of the box. The additional boundary condition, that 
<p q (z) should vanish outside the box, has to be taken into 
account also for the bound states to make them properly 
orthogonal to the unbound states. 



The wave functions tfq(z) can be found by standard 
methods^ and the energy eigenvalues are given by 



H | qk) = E qk | qk) , E qk = E q + E k 



(6) 



where E q are the energies of the QW levels. In the sim- 
plest case the QW subbands are parabolic and E k — 
h 2 k 2 /2m™ . In principle one can however introduce here 
the realistic band dispersion, although it is necessary that 
the energy function E k depends only on the magnitude 
k of k. We stress that the discontinuity in the parallel 
effective mass m\\ (see above) is a separate issue; any non- 
parabolicity in the dispersion of the basis states is due to 
coupling to the valence bands and/or higher conduction 
bands. 

Expanding in the complete orthonormal set | qk) , the 
total wave function can now be written as 



*(r)=£ J dkC q (k)\qk). 



(7) 



An important benefit of the chosen basis is that ^(r) 
fulfills the QW boundary conditions by construction. 

The spherical symmetry of the Coulomb potential is 
broken by the presence of the quantum well, and instead 
the total system adopts a cylindrical symmetry around 
the QW growth axis z. The total angular momentum 
L 2 is hence not conserved, but its projection L z on the 
growth axis is. The eigenf unctions of L z are e~ lm6 with 
m = 0, ±1, ±2, . . ., and they form a complete set. The 
same holds in Fourier space, and hence we can write 



c q (k) = c q (k,e k )= J2 



C qm (k). (8) 
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By inserting this into the expansion Eq. @ and perform- 
ing the angular integral, the eigenstates of the Hamilto- 
nian H given by Eq. Q can be expressed as 



,-im* ^ <p q (z) / kdk C qm (k) J m (kp) 



(9) 

where J m is the m-th order Bessel function. 

In principle one can diagonalize the total Hamiltonian 
H using this form of the wave function, but it leads 
to several integrals which can not be evaluated analyti- 
cally. Instead, we insert the expansion |7|. taking into 
account the angular separation according to Eq. © , into 
the Schrodinger equation H^> (r) = E^(r). As expected 
from the cylindrical symmetry, the subspaces belonging 
to different projections m are not coupled to each other. 
We can therefore solve the problem - i.e. determine the 
expansion coefficients C qm (k) - separately in each sub- 
space (for a fixed m). 

The basis expansion turns the Schrodinger equation 
into a Fredholm integral equation of the second kind. It is 
customary to symmetrize such equations by multiplying 
both sides by \fk and defining 



D qm {k) = ~\/k C qm (k) 



(10) 
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the integral equation takes the final form 



r°° r 

J o dk 1 h™,(k,k') D q>m (k') = (E — E qk ) D qm (k) - [M qql (k)+£ qq , 



D q 'rn{k^). 



where the Hermitian kernel 

h" 1 , (fc, k ) 



-e 2 Vkk> 
27re 



dz ip*(z) ip q '(z) / d9 
Jo 



-imB e -\z-z \a kk , (6) 
Ctkk' (6) 



r 



(11) 



(12) 



was calculated using the 2D Fourier transform of the 
Coulomb potential. 6 is the angle between k and k', 

and akk' (0) — V k 2 + k' 2 — 2kk' cos 9. The quantity 
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mu{z) 



- 1 



(13) 

takes care of the discontinuity of the parallel mass m\\(z), 
and the electric field £ enters through 



S, 



e£ 



dz ip*(z) z cp q >(z). 



(14) 



To solve the integral equation we discretize the fc-axis 
and approximate the integral over k by a discrete sum. 
This turns the integral equation l|ll|l into a matrix eigen- 
value problem Ti. qq , T> — ET>, where I? is a column vector 



of the (renormalized; see below) coefficients D qm {k). In 
the simplest (but most convenient) case we choose an 
equal stepsize Afc for the discretization, and the matrix 
elements are then given by 



Hf q , = (E qk 8 q<q , + M qq >{k)+£ q 
+ h™,{k,k')Ak. 



qq' J Ok.k' 



(15) 



Due to the long range nature of the Coulomb potential, 
a singularity appears in the kernel Eq. I|12l) for scattering 
in the forward direction. We do not explicitly consider 
screening in this work; instead the divergence is handled 
by averaging over small scattering angles^ In the case 
k = k' = the entire kernel vanishes exactly, but for 
other values of k — k' we obtain 



oo / .27r-Afe/2fc -imB - \z-z 1 a kk (8) 2 

h™,(k,k)*—^ I dz <p*(z) Vq .{z) d6 — + (k?0) (16) 

27re q JAk/2k a k k(0) eyfir 

A word on the normalization of the wave functions is appropriate. The matrix eigenvalue problem is solved by 
standard numerical methods, which return the eigenvectors X>W anc j eigenvalues for each eigenstate i. Typically 
such eigenvectors are normalized by ^2 qk [D q m{k)\ 2 = 1. Using the real-space density function 

pW(r)= fdv 1 [^(r'W S(r-r') ^>(r) 

dk Vk j m (k P ) / dk' j m (k' P ) J2l D &(W* D v'm(k') M*)> 



I 

the normalization condition of the corresponding wave Thus the wave function tpm will be properly normalized 
function ipm(r) becomes (after discretization J dk — > if the physical coefficients D qm (k) are taken as 



E fe Afc) 

1 = J pW(r)dr = 27rA^ |4i(fc)| 2 . (18) 



U« (*) = (19) 



q k For the particular basis and QW potential profile we 
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have chosen, most integrals which appear in the expres- 
sions above can be evaluated analytically, also when an 
electric field is included, leaving only a numerical integral 
over the angle 9. The expressions are however much too 
lengthy to include here. 



III. CHARACTERIZATION OF THE 
EIGENSTATES 

The diagonalization of the matrix problem provides the 
energies of all eigenstates of the Hamiltonian Q , and us- 
ing Eq. the corresponding eigenvectors allow us to 
evaluate the wave functions and hence matrix elements 
such as optical dipole- interaction strengths. It is however 
not immediately obvious how to identify the individual 
eigenstates as localized impurity states, QW band states 
or hybridized states, and we shall now address this ques- 
tion. 

We will adopt the notation of Ref. y and denote the 
case m = by £ and all other values by II. If the impu- 
rity is placed in the center of the well, and there is no elec- 
tric field present, the system is invariant under reflections 
z — > — z, and each eigenstate will also possess a quan- 
tum label g and u for even and odd parity, respectively. 
If this symmetry is broken, the only good quantum labels 
will be to and the energy E. It is however almost always 
possible to trace the eigenstates back to the symmetric 
situation, and we will therefore use eigenstate labels such 
as Tig even when asymmetry is present. 

Our interest is particularly focused on the lowest anti- 
symmetric impurity state, and we distinguish it from the 
other states of the same symmetry by referring to it as 
the S* state. This state has been shown to be attached 
to the second QW subband, and for narrow well widths it 
becomes resonant with the continuum of the first (lowest) 
subbandi^ It was further demonstrated that the to = 
±1 or II U states do not become resonant, but instead 
remain attached to the lowest subband. We will therefore 
from now on specialize to the case to = and often omit 
the label to. 

In the absence of an electric field, we can always iden- 
tify the states corresponding to the original QW lev- 
els (the subband bottoms), from the fact that these are 
the only eigenstates with non-zero contribution from the 
k = basis states. To facilitate the further classification, 
we define for each eigenstate i, with energy Ei, a quantity 



/>oo 

d q (E l )=2 7 r dk\D^(k)\ 2 , (20) 
Jo 

which measures the relative contribution to that eigen- 
state from the g-th basis state <p q (z); hence J2 q d q — 1 
for all states. The lowest QW level is taken as q = 0. 

Studying the values of d q (Ei) for a particular eigen- 
state i now allows us to distinguish between three types 
of states: 



• Localized impurity states only have contributions 
from higher QW levels, i.e. non-zero d q only for 
levels q with energy E q > Ei. 

• Continuous QW subband states only have contri- 
butions from lower levels, E q < Ei. (This is valid 
also in the case of a discontinuous parallel effective 
mass.) 

• Hybridized (coupled) resonant states have contribu- 
tions from both higher and lower levels. 

Our method allows us to place the impurity anywhere 
in the system. Figure ^ summarizes the qualitative be- 
havior of localized and resonant impurity states, as we 
change the impurity position from the middle of the 
quantum well to a remote location in the barrier. Some 
properties displayed schematically in this figure will be 
considered in more quantitative details in Section I VII 
whereas the remainder of this section will focus on the 
properties which characterize the various impurity states. 

Any state appearing in the energy region below the 
lowest QW level must be localized. The lowest one will 
be the impurity ground state, which is split due to the 
central-cell effect (c.f. Section ITvTl . There are several T, g 
states attached to the lowest subband, and they form 
what resembles a Rydberg series, with decreasing bind- 
ing energies converging towards the lowest subband edge. 
Actually, such series of localized states appear below each 
subband, i.e. each QW level has a set of impurity states 
attached to it (c.f. Figure The binding energy of a 
localized state is therefore to be understood as the small- 
est energy required to place an initially localized electron 
into the corresponding QW subband. It is noteworthy 
that these considerations also apply to the unbound QW 
levels. Hence if the well is so narrow that there exists 
no second bound QW level, all the S u states are still 



3D continuum 

T 



Conduction- 




Shallow impurity 
ground and 
excited states in 
the barrier 



Impurity states 
are attaclied to 
each QW level 



Localized impurity 
ground state in the well 



FIG. 1: Schematic picture, presenting the general behavior 
of how the various impurity states evolve as we change the 
impurity position, for a donor placed inside or outside a het- 
erostructure quantum well. The shaded areas represent the 
energy width of the resonant states. 
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FIG. 2: The contribution di(E) of the second QW subband 
to the eigenstates in the energy range from the first QW level 
(at 10 meV) to the second one (at 40 meV). The width of 
the Si/Sio.gGeo.i well is a=50 A, and the three curves corre- 
spond to different impurity positions zq: zo/ffl=0.05 (solid), 
0.25 (dotted), and 0.45 (dash-dot); z = is the middle of 
the well. In addition to the resonances at 15-20 meV, similar 
ones resembling a Rydberg series, appear close to the bot- 
tom of the second subband (for clarity this is only shown for 
Va=0.05). 

well-defined, but attached to the three-dimensional (3D) 
continuum. 

If the quantum well is wide enough, so that the binding 
energy of the lowest resonant state X* is larger than the 
separation between the first (q = 0) and second (q = 1) 
QW levels, this state appears below the lowest QW level, 
and is therefore localized. It is, however, still attached to 
the q = 1 level. As the well width is decreased, the X* 
state will therefore eventually appear in the continuum of 
the lowest QW subband, and we obtain a Fano resonant 
state - a localized state degenerate with a continuum. 
Similar resonant states can be formed from any excited 
impurity level (except those attached to the lowest sub- 
band), for suitable well parameters. We shall limit most 
of the discussion to the X* state, but our approach can 
equally well be applied to any of the resonant states, of 
either parity. 

As long as the impurity is placed in the center of the 
well, no hybridization or coupling can take place between 
impurity states which are resonant with subbands of op- 
posite parity. The localized X* state is attached to the 
q = 1 QW level, and therefore do(E^* ) = in the sym- 
metric situation. Moreover, without any hybridization, 
d q = S q fi for the continuous states belonging to the low- 
est QW level. Coupling will however be present as soon 
as any asymmetry which breaks the parity conservation 
is introduced - such as shifting the impurity position or 
applying an electric field. The localized state is then 
"diluted" throughout a band of actual stationary states, 
whose profile is represented by a Lorentzian resonance 
curved This hybridization mixes a certain amount of 



continuum into the X* state (giving do(E^) 7^ 0), and 
the continuous states acquire a partly localized character 
(di 0). Thus the degree of hybridization is measured 
by the coefficients di(Ei), and the resonance profile is ex- 
actly given by these coefficients when plotted against the 
eigenenergies Ei in the region between the first and sec- 
ond QW levels (c.f. Figure [2J). If the electron is initially 
placed in the coupled resonant impurity state, it will au- 
toionize with a mean life-time r = H/T determined by 
the energy width T of the resonant stated We can cal- 
culate this width by fitting a Lorentzian to the resonance 
profile. 

In the barrier, the binding energies - which are still 
measured relative the corresponding QW levels - de- 
crease towards zero, although not very rapidly. As the 
impurity is moved farther away from the well, each im- 
purity state remains bound to its initial QW level, with 
vanishing binding energy. Hence, even at very large im- 
purity distances, the impurity "ground state" appears to 
be a deep state with an energy defined by the lowest 
QW level. This is counter-intuitive; we expect the im- 
purity ground state to be a shallow state bound to the 
3D continuum. However, the electron density of the deep 
impurity states is actually localized close to the center of 
the well, and not on the impurity. When the distance 
to the well increases, their wave functions become more 
and more symmetric around the center of the well, and 
assumes the shape of the QW basis functions <p q (z). At 
the same time these wave functions are extended in the 
radial direction in the QW plane, due to the smaller bind- 
ing energy. This behavior can be understood from the 
fact that the electronic density is determined by the ef- 
fective potential in the well, and when the impurity is far 
away, the tail of the Coulomb potential is much weaker 
than the QW confinement potential. At large impurity 
distances, the Coulomb potential and the well are almost 
(but not completely) decoupled, which means that the 
deep donor states essentially derive from a single QW 
level, which can be verified by studying the coefficients 
d q . Similar observations have been made by Stopa and 
DasSarma— who also use a non-variational technique to 
treat donors in GaAs quantum wells. 

As was stated above, there are also QW impurity states 
attached to the 3D continuum (c.f. the upper right part 
of Figure When the impurity is in the barrier, these 
states are the familiar ground and excited shallow donor 
states, with the electron density localized at the impurity 
position. We can now study what happens to an initially 
decoupled shallow donor, as it is moved closer to the 
quantum well. Due to the induced coupling, Rydberg- 
like series of impurity states - localized in the well - are 
pushed down from each continuum edge, i.e. from each 
QW level. The lowered symmetry due to the presence of 
the quantum well causes mixing between the hydrogenic 
n and £ levels (actually only between levels for which 
A£ is an even number—), and only m remains as a good 
quantum label. The original shallow donor state broad- 
ens due to the coupling^ but remains bound to the 3D 



7 



continuum also when the impurity finally is located in 
the well. If the electron is initially placed in the shallow 
donor state outside a quantum well, it will after some 
time tunnel into the QW. The tunnelling time is essen- 
tially given by the energy width of the resonant shallow 
state, which can be calculated within our basis expansion 
approach, or as in Ref. l36t 



IV. STRAIN AND CENTRAL-CELL EFFECTS 

In order to compare calculated results for the ground 
state binding energy of donors in Si quantum wells with 
experimental values, it is essential to take into account 
the effects of both strain and the central-cell shift. Al- 
though there would, in principle, be some central-cell 
splitting also of the excited states, this effect is expected 
to be very small4& This is especially true for the odd par- 
ity states, since they have a vanishing envelope function 
at the impurity position. Thus the energy of the reso- 
nant S* state, which in the bulk limit corresponds to the 
2p hydrogenic state, will not be affected, and we shall in 
this section only be concerned with the impurity ground 
state. 

The conduction band minimum in bulk Si lies at about 
Ai = (0,0,0.85) in units where the X point is (0,0,1), 
with six equivalent valleys. Hence the donor Is ground 
state would be sixfold degenerate. Experimentally, how- 
ever, one instead observes three levels, with binding en- 
ergies Eb + A (the non-degenerate true ground state), 
Eb + A — Ae (twofold degenerate) and Eb + A — At 
(threefold degenerate). Here Eb is the effective-mass 
binding energy, which according to Kohn and Luttinger^, 
is about 29 meV in Si. The energies Ao, Ae and At are 
positive and, in contrast to Eb, depend on the partic- 
ular impurity species. This additional contribution to 
the binding energy is known as the central-cell effect or 
chemical shift, and is in general determined from compar- 
isons with experimental values. In the case of P donors 
in Si, Ao ~ 16.6 meV (or 14.3 meV if one uses the 
effective-mass value .©6=31.27 meV from the more elab- 
orate calculations by Faulkner—), Ae ~ 13.0 meV and 
A T ps 11.7 meVia 

The chemical shift is usually considered as arising from 
a strong potential acting only very near the impurity cen- 
ter. Although, as discussed in the Introduction, a proper 
treatment is not strictly possible in effective-mass theory, 
one could try to incorporate the central-cell effect by em- 
ploying a short-range pseudo-potential. This, however, 
fails when applied to a basis expansion method, since 
this potential now interacts not with the total wave func- 
tion (as is the case in a variational calculation) but with 
the individual basis states. If one adjusts the pseudo- 
potential parameters to produce agreement with exper- 
imental binding energies for a certain well width (such 
as a very wide one, close to bulk), the parameters will 
not give any sensible results when the basis is different, 
say for a narrow well. Instead, we will here show how to 



relate the QW central-cell shift to the bulk shifts, purely 
by symmetry considerations. 

In the effective-mass approximation, the wave func- 
tions are written as <&(r) = 0(r)u(r), where the envelope 
function 4> varies slowly over inter-atomic distances, and 
the Bloch function u is rapidly varying and describes the 
region close to the atomic nuclei. In Si there are six 
equivalent conduction band minima or valleys I, which 
we index by I £ {x,x,y,y, z,~z} (x is shorthand for —x). 
Hence one must generally consider linear combinations 
u(r) — J2e a e u ei r ) °f the Bloch functions ug of each val- 
ley bottom. If we denote the (unspecified) central-cell 
potential by V, the non-zero matrix elements between 
the valley bottom Bloch functions in Si are^S 



(ue I V | ut<) 



V , 

v 9 , 

[Vf, 



I' 



(21) 



By evaluating the matrix elements (^\ s | V | 3>i s ) for the 
six Is donor states <i>^ s (r), we can relate the parameters 
Vq, Vf and V g to the observed impurity binding energies 
as follows: 



V Q |<Mr )| 2 - (3A T + 2A E - 6A )/6, 
V g |<Mr )| 2 = (2A B - 3A T )/6, 
V f |<Mr )| 2 = -A B /6, 



(22) 



where ro is the impurity position and <f>\ s is the Is enve- 
lope function. For further details, see Appendix lAl 

We now consider donors placed in a Si quantum 
well, grown pseudomorphically between (unstrained) 
Sii-^Gca; barriers. This case differs compared to bulk 
Si in two ways. First, the Si region will be strained due 
to the lattice constant mismatch between Si and Ge. Sec- 
ond, the envelope function is different, due to the addi- 
tional confinement. On the other hand, it is reasonable 
to assume that the matrix elements Vq, V g and Vf are 
unchanged, since the central-cell potential V is apprecia- 
ble only in a region much smaller than the well width. 
Hence, we can still use the values of these parameters, as 
obtained from Eq. Il22[| using the bulk binding energies 
(determined by Ao, Ae and At), but replace (f>is by the 
QW envelope function. 

By taking into account the effects of strain and the 
central-cell shifts simultaneously, one finds - in addition 
to an overall shift - that for (OOl)-grown wells, the six- 
fold degenerate donor ground state splits into five energy 
levels, whereof one is twofold degenerate— The details 
of this calculation are presented in Appendix \B\ The lat- 
tice constant of unstrained Sii^Ge^ bulk alloy has been 
parametrized^ from experimental dataS as 

a(x) = o(Si) + 0.0200326x(l - x) + [a(Ge) - a(Si)]a; 2 , 

where < x < 1 and a(Si) and a(Ge) are the lattice con- 
stants of bulk Si and Ge (in nanometers) . Inserting this 
into Eqs. 1|B4J) and (|B5J| we find that the strain compo- 
nent parallel to the QW interfaces en > 0, and hence the 
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strain parameter £ = — S u e|| (2C12/C11 + l)/3 denned in 
Eq. IBlOj l is negative. Here S u is the deformation poten- 
tial and Cy are the components of the stiffness tensor. 
Even with only a small content x of Ge, the condition 
that the strain splitting is large compared to the central- 
cell shifts (£ » A ) is fulfilled^ Therefore, from the 
results in Appendix [5] the lowest valleys will be Aj_ , 
perpendicular to the QW interfaces, to which two donor 
states are attached. Four other donor states are associ- 
ated with the higher A|| valleys, which are strain split 
from Aj_ by 3|£|. 

In our calculations we assume that the strain splitting 
of the conduction band is included in the band offset U, 
and consider only the two lowest donor states, belong- 
ing to Aj_. These two states are shifted from the QW 
effective-mass value Ef, by 



|^Qw(rp)| ; 
|<Mr )| 2 



-An 



2A E /3 
A T 



(23) 



The Bloch functions of the two states are the symmet- 
ric and anti-symmetric combinations of the two A^ val- 
leys. The quantum well envelope function amplitude 
|'0Qw( r o)| 2 is obtained from the basis expansion coef- 
ficients using Eq. I|17|l . We will return to the question 
how to correctly determine the bulk amplitude |0i s (r)| 2 
in Section IVll 

The central-cell shifts are negative, which means that 
the binding energies are increased. The impurity ground 
state now corresponds to the (non-degenerate) 2 As/3 
state, since Ae — At for typical donors (P, As, Sb) 
in Si4£ It is noteworthy that, since furthermore A^ — 
Ao, the central-cell effect in the quantum well case can 
be much smaller than in bulk Si, at least if the ratio 
of the wave functions is not too large. In bulk Si, the 
donor ground state is shifted by Ao = 17 meV, but in 
the QW the ground state shift for P donors is only Ao — 
2Ae/3 w 8 meV, assuming that the envelope function 
ratio in Eq. (|23|l is of the order unity. On the other 
hand, this ratio is expected to be above unity for narrow 
wells, due to the additional confinement, and so in this 
case the central-cell effect further increases the already 
enhanced ground state binding energy. 



V. OPTICAL ABSORPTION 

If the impurity is placed exactly in the middle of the 
well (and no electric field is present), the selection rules 
prohibit the radiative decay from the lowest subband to 
the impurity ground state by a dipole transition with po- 
larization parallel to the QW growth direction. However, 
if some asymmetry is introduced, the transition becomes 
allowed from the part of the subband which is hybridized 
with the resonant anti-symmetric E* impurity state. 

In the dipole approximation, Fermi's Golden Rule gives 
the probability per unit time of optical absorption, at the 



frequency u>, between states i (initial) and j (final) as 
27r / eA 



Wji = -r — |0|e. P | ? )r^-^-M (24) 
n \ mc J 

where Aq is the magnitude of the electromagnetic vector 
potential, e the photon polarization vector, p the mo- 
mentum operator, and Ei the energy of eigenstate i. 

Assuming that the system contains N independent im- 
purities, we may relate the absorption rate Wji to the 
absorption coefficient A as 



A(huj) = 



N 2-irhc 

n klo\Aq 



2 A^i 

3 



311 



(25) 



where O is the total sample volume and n the refraction 
index. In a quantum well, a more natural quantity to 
consider is the absorption cross-section 



An 

~N~ 



4n 2 h 2 c 



E 

3 





e ■ p 












m 





6(hujji — Hu) 



where 



Ei)/h and a = e 2 /hc is the 



(26) 
fine- 



U)ji = (Ej 
structure constant. 

The summation in Eq. (j26|) in principle involves, 
among other things, final states j in all valleys with any 
spin polarization. We can however ignore the four A|| val- 
leys altogether, since the strain splitting is much larger 
than the impurity binding energy. We shall consider ab- 
sorption only from the ground state, and so the initial 
state i is an equal mix of the two Aj_ valleys, as dis- 
cussed in Section llVl with fixed (but arbitrary) spin. The 
dipole operator e ■ p does not affect spin, and connects 
only the parts of the initial and final states belonging to 
the same valley. This holds also when the temperature is 
large enough that the occupation probabilities of the two 
states 112 and u§ (c.f. Appendix |B|) are roughly equal; in 
this case we must also average over the two states. Thus 
the sum over j can be restricted to a sum over the final 
energy E and the different cylindrical subspaces. It is 
furthermore enough to evaluate the matrix element be- 
tween the envelope functions only, as long as we replace 
the electron mass m by the effective mass m* in the di- 
rection parallel to p4£ 

Since one would most naturally place the donors inside 
the QW for optical applications of this kind, we assume 
that the relevant effective masses are those of the well, 
and write 



e • p cos 6p x 



(27) 



where 9 is the angle between the photon k- vector and the 
QW growth direction z. The light is assumed to be plane 
polarized in the plane spanned by the normal to the QW 
plane and the photon k-vector. The ir-axis is defined 
as the intersection of this plane and the QW interface 
planes. 
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In principle Eq. I|27(l allows for interference between 
the two terms, when inserted into Eq. J2HJ). The in- 
terference terms are essentially (i \ p x \ j) (j \pz\i) or the 
complex conjugate thereof, to be summed over all final 
states j. Here i is the impurity ground state, which 
belongs to the m, = cylindrical subspace. The sum 
over final states can be split into two parts, where we 
first consider all states j in a particular subspace rrij, 
and then sum over all subspaces. The first factor is 



\l\Px \J) oc mj , mi ±i 



which follows from the selection 



rule Am = ±1 (see below) for the p x operator. On the 



other hand, (j \p z 



J rrtj ,771; 



since Am — for p z 



Thus, the interference terms vanish, and we may write 

cr = a x cos 2 9 + a z sin 2 9, (28) 

which defines a x and a z . As we will see in Section ED 
the absence of interference is a contributing factor to pro- 
ducing a symmetric absorption peak, contrary to what is 
often observed for optical transitions involving resonant 
states32i 

Let us first consider z-polarization. To eliminate the 5- 
symbol in l|26|) we turn the sum over final states j into an 
integral. However, we have to consider that the relevant 
density of final states in this case is not the usual 2D 
density of states. If we double the "size" of our system 



(Afc — > Afc/2) we only get twice as many eigenstates, 
since the discretization in Section [H] is carried out over 
the magnitude fc of k. Hence the density of states in 
fc-space is dn/dk = 1/Afc, and in energy space 



1 



dn dn dk 



dE k dk dE k hAk V 2E k 



(29) 



where the parabolic dispersion relation E k = h k / 2m™ 
was again assumed. The presence of Afc will properly 
renormalize the matrix element when this is evaluated 
from the discretized expansion coefficients (c.f. the dis- 
cussion on normalization in Section[H|. Now we may per- 
form the integral over energy to remove the <5-function, 
and the absorption cross-section can be evaluated at 
u) = ujji for all eigenstates j. The result is 



27r 2 a 



3/2 



(j \Pz 



(30) 



By using the form Eq. 10 of the wave function for the 
states i and j (belonging to the subspaces m and m', 
respectively), the remaining matrix element becomes 



J 



(j \Pz 



271 



r 



qq' 



(31) 



Since the integral yields 2irS mtm i, we obtain the selection 
rule Am = for this polarization. The integral over fc can 
be evaluated from the matrix eigenvectors U, taking the 
normalization relation Eq. (|19|l into account. We can fur- 
thermore use the commutator [Hqw,z] = —ihp z /mj_(z) 
and the fact that the states q and q 1 are eigenstates of 
Hqw, defined in Eq. J2J), to rewrite 

i f°° 

W \Pz\q) = T\ E q' - E q) I ¥* q ,{z) Z m ± (z) tf q (z) dz. 
" J — oo 

Next we consider x-polarization, for which the selec- 
tion rule is Am = ±1, as will be demonstrated shortly. 
We will focus on the energy region between the first and 
second QW subbands, where the lowest resonant state S* 
appears. Since there are no resonant states with m = ±1 
in this region, and since the ground state belongs to the 
m = subspace, we can ignore mixing of the QW sub- 
bands and represent the final states as normalized plane 
waves belonging to the first subband q = 0: 

. . exp(ikW ■ p) 

b) = z ^ 

where C is the normalization length. 



)(*), 



(33) 



To derive the selection rule, one may rewrite the mo- 
mentum operator p x in cylindrical coordinates and act 
on the cylindrical expansion of a plane wave 



J m (kp), where $ is the angle between the 



vectors k and p. One then obtains 



( j I Px 



nhk x 3) 
C 



oo 

E 

rn=—o 



3m. ±1 



(34) 



which shows the selection rule explicitly; the summation 
is trivial and gives a factor of 2. As before, i refers to the 
donor ground state, represented by an expansion of the 
form @ . It is not surprising that the momentum opera- 
tor p x picks up the cc-component of the final wave vector 
k") ; the result Eq. (|34J) can of course also be obtained by 
acting with p x on the plane wave e lk p without changing 
to cylindrical coordinates. But in that case the selec- 
tion rule will not appear, since it is not actually present 
in the final result, after the summation in Eq. Q34JI has 
been carried out. 

We remove the 5-function in Eq. I|2()|) by integrating 
with the usual 2D density of states per spin £ 2 m^ /2irh 2 . 
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In result, we arrive at 

27r 2 q HKq=oj (») , , 2 , 

° x = ^Kk^- P Uo(^=o„) (35) 

where K q j is defined from the energy conservation rela- 
tion 

Ei + hua = E q + -^-, (36) 

where Ei is the energy of the donor ground state and E q 
the energy of the q-th. QW level. The expansion coeffi- 
cient V q l \K q j) can be found by interpolating the matrix 
eigenvectors over k. 

VI. RESULTS 



huge (typically N will be between 1000 and 2000). The 
size of the matrix problem which can be diagonalized is 
a limiting factor, and it is not always possible to obtain 
optimal accuracy, especially for very wide wells. Never- 
theless, we have made sure that the errors in the results 
presented here are less than 1%. The only exceptions are 
the energy width of the resonant state when the impurity 
is deep inside the barrier (in which case the width itself 
is very small), and the ground state wave function am- 
plitude; here the errors may be up to 5%. However, the 
accuracy in the ratio of the wave function amplitude to 
the corresponding bulk value - which really is the quan- 
tity of interest - is much better than the accuracy of the 
amplitude itself. The wave function amplitudes are also 
more sensitive to the size of the outer box C z (which nor- 
malizes the " continuum" basis functions) than the eigen- 
values are. Again, the convergence can be controlled by 
increasing this size as much as needed. 



A. Numerical aspects 

In this section we present numerical results for shallow 
Coulombic donors in (OOl)-grown Si/Sii_ x Ge x quantum 
wells. The electronic parameters (effective masses, band 
offsets and deformation potentials) of these systems are 
known from the systematic study by Rieger and Vogli^i 
We make the approximation, as discussed in Section ITU 
that the difference in dielectric constant can be ignored, 
and use that of bulk Si throughout. Furthermore, al- 
though the non-parabolicity of the subbands in strained 
Si quantum wells has been found to be considerable, we 
take the basis subbands to be parabolic, since the non- 
parabolicity parameter is not well known^ For small 
contents x of Ge the effective masses are very similar for 
the well and barrier regions, and therefore the disconti- 
nuity in the parallel effective mass turns out to have no 
observable influence on the results for these particular 
systems. 

For the matrix problem to be of finite order we must 
limit the number of QW levels to include in the basis, 
and also the integration (sum) over k must be cut off 
somewhere. This is the only real approximation in the 
method (other than those inherent in effective-mass the- 
ory). Therefore, in order to assure the numerical ac- 
curacy of our results, we always include all bound QW 
levels, and increase the number of unbound levels un- 
til the eigenvalues do not change. The same procedure is 
applied to determine the required density and range of k- 
values. The contribution from the unbound states drops 
off rapidly with increasing energy, although for narrower 
wells it is naturally important to include more unbound 
states. 

Quite a large range of fc-values is needed to reach con- 
vergence, especially in the ground state wave function 
amplitude, since the expansion coefficients decay rather 
slowly, roughly as k~ 3 5 . This means that, since at the 
same time one wishes to keep Afc as small as possible 
for accuracy, the N x N matrix problem becomes fairly 



B. Binding energy and resonant state energy width 

In agreement with variational calculations, we find that 
the binding energies are generally increased due to the 
additional confinement presented by the quantum well 
potential, both for the ground and the excited states (c.f. 
Figures |21 and 0J. Furthermore, since the E* state is 
attached to the second QW level, it appears at higher 
and higher energies as the level separation increases with 
smaller well widths. This holds until the well becomes 
very narrow, when the binding energies decrease again, as 
shown in Figure [21 When the well width is much smaller 
than the radius of the impurity state, one can view the 



60 




10 -J 



1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 

50 100 150 200 250 300 

Well width (A) 

FIG. 3: The binding energy of the ground state (solid curve) 
and the state (dashed curve) for a Si/Sio.sGeo.2 quantum 
well of varying width. The impurity is located in the middle 
of the well, and the central-cell shift is not included. For wide 
wells the binding energies approach the bulk values (29 meV 
and 7 meV, respectively), indicated by the dotted lines. The 
inset shows the energy separation between the S* state and 
the ground state (GS). 
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FIG. 4: The impurity ground state binding energy (left-hand 
vertical axis) for Si/Sio.sGeo.2 wells of widths (from top to 
bottom) a=30 A, 40 A, 50 A and 60 A, and for the same 
widths and symbols the position of the resonant state 
(right-hand vertical axis). The right-hand vertical scale is 
normalized so that corresponds to the first QW level, and 1 
to the second one. The QW level energies are, from narrowest 
(30 A) to widest (60 A) well: £ =26.8, 17.1, 11.8, 8.7 meV 
and _Bi = 102.7, 67.1, 46.8, 34.4 meV. The vertical line marks 
the edge of the well, and the impurity offset is measured from 
the middle of the well. The central-cell shift is not included. 



system as an impurity situated in the 3D bulk, slightly 
perturbed by the narrow quantum well. For a vanishingly 
narrow well the binding energies therefore approach the 
values corresponding to the bulk barrier material. More- 
over, below a certain well width there is no bound odd- 
parity QW level, and the E* state then becomes attached 
to the continuum, which does not move as the well width 
is changed. These effects on the ground and resonant 
state binding energies have also been observed in varia- 
tional calculationsiSiii 

For wide wells, the ground state binding energy con- 
verges to the effective-mass value of the bulk hydrogenic 
Is donor state in Si at a width of about 250 A, which 
is faster than for GaAs/Al^Gai-^As quantum wells. 9 
Clearly, the reason is the smaller radius of the impurity 
states in Si compared to GaAs. By the same reasoning, 
the convergence of the excited E* state to the correspond- 
ing bulk state 2po is slower than for the ground state, 
since the smaller binding energy of the excited state cor- 
responds to a larger radius of the state. From Figure 01 
we further note that for uniformly doped wells, one effec- 
tively obtains broad bands of impurity energies, instead 
of a set of sharp levels. The band-width of the ground 
state can be almost comparable to the spacing of the QW 
levels, in particular if the central-cell effect is also taken 
into account (see below). 

The decay of the binding energies in the barrier is very 
slow. It follows a power-law dependence, but there does 
not appear to exist any common exponent; for the curves 
shown in Figure 0] the exponents are approximately 0.7- 
0.8 for the ground state and 0.2-0.4 for the resonant 
state. Nevertheless, at very large distances the impurity 
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FIG. 5: The energy width Y of the resonant S* state for the 
same quantum wells as in Fig. 21 with the same symbols for 
the curves. 



states coalesce with the quantum well levels; see the dis- 
cussion in Section IlIII regarding the impurity wave func- 
tions when the donor is placed in the barrier. 

Figure |S] shows the energy width T of the resonant E* 
state. As the impurity is moved away from the center of 
the well, the width increases due to the enhanced cou- 
pling (the increased asymmetry); this could also be seen 
from Figure When we continue to move towards the 
barrier, T reaches a maximum value at about 35% offset. 
After that it decreases, due to the reduced overlap be- 
tween the two lowest QW basis functions. In the barrier, 
the width continues to decay according to a power-law 
dependence with an exponent (which is not common for 
different well widths) of the order 2-3. This is in con- 
trast to the resonant state formed from the shallow donor 
states when the impurity is in the barrier (shown in the 
upper-right part of Figure for which the width decays 
exponentially^ 

For the same relative impurity offset, the resonance 
is broader in wider wells, which is an effect of the fact 
that the E* state appears closer to the bottom of the 
band. In Figure|2]the same can be observed for the higher 
excited E„ resonant states, which are narrower the higher 
energy they have. We furthermore see from Figure El 
that a widening of the resonance causes its amplitude 
to decrease. The smaller the energy width is, the more 
localized the impurity state is, but the effect on scattering 
and optical properties can still be pronounced since the 
amplitude of the resonance at the same time is larger. 

The behavior of the resonance position for the widest 
well (a = 60 A) is slightly different from the other curves 
in Figure 0] For this well width, the resonance is close 
to the subband bottom for central impurity positions. 
Therefore, the resonance profile becomes rather asym- 
metric, since it cannot be continued below the lowest 
subband. The amplitude of the resonance, as measured 
by d\ ) , is furthermore very small as long as the res- 
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onance is close to the band bottom. At the same time 
the resonance width T is the largest for this well width, 
at intermediate impurity offsets. To accommodate this 
width, the resonance position shifts away from the sub- 
band bottom more quickly in this case than for the other 
well widths. 



C. The central-cell effect 

The results discussed above for the ground state bind- 
ing energy indirectly contain the effect of strain, since this 
gives rise to the QW band offset, but do not include the 
central-cell effect, and therefore represent the effective- 
mass binding energy. To take the central cell into account 
for the two donor states attached to the lowest A_l valleys 
is however straightforward, by using Eq. (J2SJ). The bulk 
shifts Ao, Ae and At depend on the particular impurity 
species (values for phosphorus are given in Section Hvjl. 
and in addition we need to know the ratio of the QW and 
bulk envelope functions at the impurity position. 

Regarding the bulk envelope function, one could imag- 
ine using a variational wave function, such as the nor- 
malized non-isotropic "hydrogenic" function 

^ (p ' z) = vhi exp (~ V?+?) ' (37) 

originally used by Kohn and Luttingen4ii to obtain the 
value 29 meV for the effective-mass binding energy in 
Si (the energy is minimized by a — 2.478 nm and b 
1 .420 nm). However, although variational functions may 
produce rather accurate estimates of the energy, they do 
not, in general, give a correct picture of the wave func- 
tion. This is illustrated by FigureEJ where we see that in 
the region close to the impurity, the variational function 
ifiis (p, z) given by Eq. 137J1 may differ substantially from 
the wave function tpQwip, z ) obtained with the basis ex- 
pansion method, although they both reproduce the same 
binding energy. The two envelope functions do however 
agree very well in the exponentially decaying tail, which 
is precisely the part which determines the binding energy. 

In the inset of Figure [7J the value of the quantum well 
envelope function amplitude, at the impurity position, is 
plotted as a function of the well width. Comparing with 
Figure |21 we see that at a width of 250 A, the ground 
state binding energy is extremely close to the variational 
value, which is reasonable since the well is much wider 
than the radius of the impurity state (given by a and 
b). However, the QW envelope function amplitude differs 
substantially from |<pi s (0,0)| 2 = 0.0365 nm~ 3 . Instead, 
as the well becomes wider, the amplitude converges to a 
value of approximately 0.023 nm~ 3 , which we therefore 
will take as the value for the bulk amplitude |</>i s (r )| 2 . 

Once the bulk amplitude thus has been determined, we 
can use the ratio of the QW and bulk envelope functions, 
shown in Figure [7| to evaluate the chemical shift for any 
wells width and impurity position. When the impurity is 
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FIG. 6: The surface plot is the difference |<pis| 2 — IV'Qwl 2 
for a wide (200 A) Si/Sio.sGeo.2 quantum well, and the inset 
shows the QW envelope function |i/>qw(.z)| 2 (solid) and the 
bulk variational function |<pi s (z)j 2 (dashed) for p = 0. Both 
functions are normalized and give the Si bulk binding energy 
(29 meV; c.f. Figure 



close to the center of the well, the QW wave function in 
narrow wells is strongly enhanced at the impurity posi- 
tion, whereas for wider wells it approaches the bulk value. 
On the other hand, moving the impurity to the edge of 
the well leads to a rapid reduction of the ratio. It was 
mentioned in Section lTTTI that when we place the donor in 
the barrier, the wave function of the deep impurity states 
is localized in the well, and not on the impurity. Now we 
find that even when the impurity is still inside the well, 
the wave function maximum does not exactly coincide 
with the impurity position for asymmetric locations. By 
studying the probability density one finds that the wave 
function is pushed towards the center of the well. One 
can view this as the wave function being reflected away 
by the barriers. 

For wells wider than 50 A, the ratio of the bulk and 
QW wave function amplitudes is close to unity, and 
hence, as was mentioned in Section llVl the central-cell 
shift is substantially smaller in the quantum well than in 
bulk Si. For off-center impurity position this applies to 
even narrower wells (c.f. Figure 0). Note also that for 
very wide wells, the shift as calculated from Eq. \\TM is 
different from that of bulk Si, due to the strain. Thus 
we can expect that the central-cell shifts are smaller in 
strained bulk Si than in the unstrained material. Since 
the chemical shift is much larger for impurities placed in 
the center of the well than for positions close to the bar- 
rier, the effective band of impurity energies mentioned 
above is furthermore widened by the central-cell effect. 
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FIG. 7: The ratio of the QW envelope function amplitude to 
the bulk value, at the impurity position ro, for three different 
Si/Sio.8Geo.2 wells of widths (from top to bottom) a=30 A, 
50 A, and 70 A. The vertical line marks the edge of the well, 
and the impurity offset is measured from the middle of the 
well. The bulk amplitude is defined as the wide-well limit 
(0.023 nm~ 3 ) of the inset plot, which shows the QW envelope 
function amplitude as a function of the well width. The am- 
plitude closely follows the variation of the ground state bind- 
ing energy (c.f. Figure [3J , also for very narrow wells when 
both the binding energy and the amplitude again approach 
the bulk value. 



D. Electric field and optical absorption 

When a static electric field is applied across the quan- 
tum well, the parity symmetry is broken even when the 
impurities are placed in the center of the well, as in the 
case shown in Figure |S1 The energy width T of the reso- 
nant E* state - and hence the degree of coupling - can 
be controlled by varying the electric field. Also the tran- 
sition energy from the resonant state to the ground state 
can be fine-tuned in the same way, but the tuning range is 
small compared to the resonance width. By varying the 
QW parameters, on the other hand, this energy can be 
tuned over a vast range, as shown in the inset of Figure[3] 

As the electric field is increased from zero, the reso- 
nance width naturally also increases due to the induced 
coupling between the resonant state and the continuum. 
Still, the detailed behavior of the curves in Figure 03 war- 
rants further comments. If we take the energy of the zero- 
field conduction band bottom as a fixed reference, we can 
study the field-dependence of the energy positions of the 
impurity states and the quantum well levels (which are 
also shifted by the presence of the field). The positions 
of the second QW level and the resonant E* state - both 
of which have odd parity without any applied field - are 
practically unaffected by the field. On the other hand, 
the lowest QW level and the impurity ground state (even 
parity states) are deflected downwards. Hence, as the 
field is applied, the resonance appears at higher energies, 
relative to the subband bottom. We noted above that 
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FIG. 8: The solid curve shows the energy width V of the 
resonant E* state, as we apply a static transverse electric 
field to a Si/Sio.sGeo.2 well of width 55 A. The impurities are 
placed in the center of the well. Also the energy separation 
(dashed curve) between the resonant state and the ground 
state (GS) is presented. 



the energy width is smaller for higher-energy resonances, 
and the interplay between this effect and the increase in 
the width from the stronger coupling gives rise to the 
plateau at about 50-70 kV/cm. At yet higher fields, the 
enhancement of the coupling will however dominate and 
the width increases again. Alternatively, one can con- 
sider that the electric field "pushes" the wave functions 
towards the side of the well. This affects odd and even 
states slightly differently, which influences the overlap 
between the first and second QW levels, and this is an 
essential factor in determining the coupling and the res- 
onance position. 

Finally we have also calculated the absorption cross- 
section a (c.f. Figures 151 and ITU)! from the donor ground 
state, in the case when the E* state is coupled to the first 
subband by placing the impurity asymmetrically in the 
well. According to the dipole selection rules, discussed 
in Section there is no absorption from the ground 
state to the unperturbed first subband (q = 0) if the 
incoming light is linearly polarized along the QW growth 
direction. However, when the E* state hybridizes with 
this subband, it mixes a certain amount of the second 
QW subband (q — 1) into the continuous states of the 
first subband, and thus <r 2 ^ 0. Hence in an energy 
region, determined by the width T, around the resonant 
state, absorption is allowed even at normal incidence (9 = 
90°). This is seen in the figures as the narrow peak at 
about 55 meV. 

The shape of the normal incidence absorption peak is 
symmetric. This follows from the fact that there is no 
interference between the matrix elements for absorption 
into the localized and continuous parts of the hybridized 
resonant state, respectively (c.f. Section 01 . Given this, 
and the energy-independent 2D density of states, the ab- 
sorption spectrum will have the same shape as the reso- 
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FIG. 9: The profile of the absorption cross-section a for a 
40 A wide Si/Sio.ssGeo.is quantum well, with the impurity 
offset 4 A from the center of the well. The electron is initially 
assumed to be in the impurity ground state; the central-cell 
shift has however not been taken into account for the photon 
energy scale. For the actual values of the cross-section, see 
Figure E3 The coordinate system for the incidence angle 9 is 
defined in Section W\ 
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FIG. 10: For the same system as in Figure^] we here show the 
absorption cross-section for some selected incidence angles 9; 
from top to bottom 9=0°, 30°, 60°, and 90°. 



nance profile, which is symmetric in the displayed (and 
nearly all other) cases, and the asymmetry often observed 
for resonant states22iM is hence absent. 

For angles 9 < 90°, absorption is also allowed for po- 
larization parallel to the QW interfaces, given by a x , and 
this gives rise to the broad background peak. Further- 
more, at about 85 meV photon energy, absorption into 
the second QW subband also becomes possible; in this 
case the selection rules allow absorption for any polar- 
ization. The cross-section is finite at the second subband 
edge, reflecting the step-like 2D density of states. 

The amplitude of the resonant state absorption peak 
is not very large compared to the background, except for 
near-normal incidence (6 ~ 90°). Nevertheless, the ab- 
sorption cross-section a z is still of the same order as the 
cross-section for impurity absorption in bulk Si^ It is 
expected that the considered intra-impurity transition is 



particularly strong, since the transition is analogous to 
the so-called resonance line (2p — > Is) in atomic hydro- 
gen. 



VII. SUMMARY AND DISCUSSION 

In this paper we have presented a unified approach for 
calculating the energy levels of shallow donors in het- 
erostructure quantum wells. By turning the Schrodinger 
equation - containing both the QW profile and the im- 
purity potential - into a matrix eigenvalue problem, we 
obtain a complete description of the entire energy spec- 
trum for all donor positions, both inside and outside the 
well. 

Applying this method to Si/Sii-^Gcj; quantum wells, 
we have calculated the binding energies of the ground 
state and the lowest anti-symmetric resonant state for 
several well widths and impurity positions. The depen- 
dencies of the binding energies on the Q W parameters fol- 
low the same general behavior found in variational calcu- 
lations of GaAs/Al^Gai-^As systems. When the donor 
is placed asymmetrically in the well, or a transverse elec- 
tric field is applied, the resonant state can hybridize with 
the continuous subbands, and we can within our method 
evaluate the resonant state energy width, which is di- 
rectly related to the life-time. 

Si/Sii-xGez quantum wells are strained due to the 
lattice mismatch, and this was investigated in detail, 
along with the central-cell effect. By expressing the 
QW central-cell shift through the parameters of the split 
donor states in bulk Si, and the amplitude ratio of the 
QW envelope function to the bulk one, it was shown that 

- depending on the impurity position and the well width 

- the QW shift can be both smaller or larger compared 
to the bulk case. 

We have also compared the commonly used variational 
function for donors in bulk Si, to the ground state wave 
function obtained with the basis expansion. Even though 
both functions give the same energy if the well is wide 
enough, they are far from identical. No assumption re- 
garding the shape of the wave function is made in our 
non-variational method, and we therefore conclude that 
the " hydrogenic" function is not capable of giving a cor- 
rect description of the donor wave function. Instead we 
obtain the correct envelope functions of the localized, res- 
onant and continuum eigenstates, and can then evaluate 
various matrix elements. As an example we present the 
ground state optical absorption spectrum, which shows a 
strong dependence on the direction of the incident photon 
compared to the quantum well axis, due to the selection 
rules. 

The possibility to populate the resonant state by elec- 
trically pumping the electrons in the QW subbands, 
could be used to create a conduction band resonant state 
laser in a Si quantum well. Provided the carriers can 
reach the resonance energy without being scattered by 
other mechanisms, they can - through the hybridization 
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- be captured into the localized part of the resonant state. 
They may then make an optical transition to the impurity 
ground state, a transition which is particularly strong, 
or to some excited localized state. Since, as was shown, 
the impurity states are attached to the QW levels, it is 
possible to tune the intra-impurity transition energies by 
varying the well parameters. The tuning range for the 
ground state transition was shown to extend from 25 to 
150 meV by changing the Si well width between 2 and 
15 nm (the central-cell effect increases this range some- 
what). Additionally, some fine-tuning could possibly be 
achieved by applying a transverse electric field. 

For the Si wells we have considered, the resonance 
width can be as large as 10 meV, which gives a very short 
life-time of about 0.1 ps. One may therefore expect that 
the resonant states can have a strong influence not only 
on the mobility, from the pronounced resonant scattering 
mechanism which appears in a narrow energy region, but 
also on the noise spectrum, due to the capture and re- 
emission process. These features are in fact present even 
when the impurities are placed outside the well, since 
there is coupling also between the shallow barrier donor 
state and the QW subbands. 



VIII. ACKNOWLEDGEMENTS 

This work was performed within the Nanometer Con- 
sortium at Lund University, and was supported by grants 
of the Swedish Foundation for Strategic Research, the 
Swedish Research Council (TFR-THZ 2000-403), NorFA 
(000384), the Russian Foundation for Basic Research, the 
Russian Academy of Science, and the Russian Ministry 
of Science. I.N.Y. is grateful for the support from the 
Wenner-Gren Foundation. 



Using the valley bottom Bloch functions ui as a basis, 
the Bloch functions Ui can be represented as vectors of 
the coefficients a\. 

However, we also know that the six Is states transform 
as the group, and hence their Bloch functions it, can 
be written as^i 

^ = (1,1,1,1,1,1)^/6, 

v$ = (l,l,-l,-l,0,0)/2, 

uf = (1,1,1,1, -2, -2)/2V5, 

U W = (1,-1,0,0,0,0)^/2, 

4 2) = (0,0,1,-1,0,0)/V2, 

uf> = (0,0,0,0,1,-1)/^. 

The label A refers to the non-degenerate ground state 
with energy — Ao relative to the effective-mass value, E 
is the doublet state with energy — A + Ae, and T is the 
triplet state with energy — Ao + At- 

Since the envelope function is assumed to vary slowly 
over distances comparable to the range of the central-cell 
potential V, we can write 

<*i.|^|$i.) = |0 ls (ro)| 2 K|^l^> (A4) 

where i"o is the position of the impurity. The overlap 
matrix elements between Bloch functions from different 
valleys I are given in Eq. I|21|) . By inserting the Bloch 
functions given in Eq. (|A3|) into Eq. I|A4(1 , we can express 
the donor energies in terms of Vb, V g and Vf as follows: 

(* A | V \ * A ) = |<Hr )| 2 (Vb + V g + 4V» = -A , 
(*b|V| * b > = l<Xro)| 2 (V + V g - 2V f ) = -A + A E , 
(*r | V | * T > = I^MI 2 (Vb - V B ) = -A + A T . 
Inverting these relationships gives the results Eq. 12 I'D . 



APPENDIX A: RELATING THE CENTRAL-CELL 
PARAMETERS TO THE BINDING ENERGIES 

In the effective-mass approximation, the wave function 
close to each conduction band minimum, i.e. in each val- 
ley £, is assumed to be of the form ^(r) = <j>(r)ue(r), 
with a common envelope function <p. Here Ui is the 
Bloch function of the valley bottom, which can be written 
Ui(r) — exp(ik^ • r)Uk e (r), where k^ is the wave vector 
of the respective conduction band minimum, and U is a 
function with the periodicity of the lattice. 

In bulk Si the six conduction band minima £ £ 
{x,x,y,y, z,z} are degenerate. Therefore we write the 
wave function of each Is impurity state as 



(rK(r) (i= 1.--6), 



(Al) 



APPENDIX B: DONOR STATES IN STRAINED 
SILICON 

The strain Hamiltonian for the conduction band in Si 
can be written^ 

Strain = S d Tr(e) + E u (k ■ e ■ k) (Bl) 

where Tr means the trace and e is the strain tensor. The 
constants and S M are the deformation potentials, and 
k is a unit vector along one of the equivalent valleys 
{x,x,y,y, z,z} in the unstrained material. In the basis 
of the valley bottom Bloch functions (c.f. Appendix lA")l . 



strain 
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The second term represents an overall shift and can be 
ignored altogether. Here 1 is the 6x6 unit matrix. 

For a pseudomorphically (OOl)-grown strained layer, 
the strain tensor takes the form 



(B3) 



where e 1 1 and e± are the strain tensor components par- 
allel and perpendicular, respectively, to the interface 
planes. These can be expressed in terms of the unstrained 
lattice constants of the layer ai and the substrate a s , 
by evaluating the new lattice constants in the strained 
layes^ 



To abbreviate the expressions we have introduced 




(B4) 



where Cy are the components of the stiffness tensor. We 
then have 



an a±_ 2Ci2 
'|| = li e ± = 1 = --=r~ e 



(B5) 



from which wc sec that the two components have opposite 
signs. 

In order to take the effects of strain and the central cell 
into account simultaneously, we make a unitary transfor- 
mation of Strain to the basis given by the states ua , Ujs 
and Uj) defined in Appendix ^ and add the central-cell 
contribution 
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which naturally is diagonal in this basis. 

The resulting Hamiltonian matrix is, apart from a com- 
mon diagonal constant — Ao, 

/ -V2£ \ 












A, 







A 7 



A, 



\ A T + 2£/ 

(B7) 

which is easily diagonalized; the once sixfold degenerate 
donor ground state splits into four non-degenerate and 
one twofold degenerate states (Ref. |4ll §37) , with ener- 
gies 



£2,3 

£4,5 
£6 



-Ao- 
-A - 

-Ao- 
-A 



A E -£, 

At — £,i 
A T + 2£. 



2x + 9x 2 



(B8) 



£ = H„(e_L - e||)/3 



(B9) 



The Bloch functions of the six eigenstates are££ 

u x = (1,1, -1,-1, 0,0)/2 
u 2 = (auA + u% )/ \/l +a 2 
u 3 = (Bi-atifij/Vl + a 2 

r- 

u 4 = (1,-1,0,0,0,0)/V2 
u 5 = (0,0,1,-1,0,0)/V2 
u 6 = (0,0,0,0,1,-1)/V2 

expressed in the valley bottom Bloch function basis (c.f. 
Appendix ^| with ua and given in Eq. i|A3(l , and 



1 + x + Vl + 2x + 9x 2 ' 



(B12) 



If the strain splitting S n (ej_ — e 1 1 ) ^ s mucn larger than 
the central-cell splitting Ae, as is the typical situation 
in Si grown on a Sii^Ge^ substrate x 3> 1 and 
a ~ — 1 / V2 ■ In this limit the eigenfunctions u 2 and 
«3 become 



u 2 « (0,0, 0,0,1, 1)/V2, 
U3« (1,1, 1,1,0, 0)/2, 



(B13) 



with energies 



e 2 « 2£ - A - 
£3 ~ -t, ~ A 



2Ab 
3 

Ab 
3 



(B14) 



Thus in this case the six states separate into two groups, 
where u 2 and uq are comprised of the two & z -valleys, and 
the four other states are not coupled to these valleys at 
all. 

In the opposite limit, if the strain is very small, i«l 
and we see that u\ and u 2 originate from the doublet and 
w 3 from the ground state. 
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